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Abstract 

Background: Drug resistance displays a problem for the therapy of Mycobacterium tuberculosis infections. For 
molecular resistance testing, it is essential to have precise knowledge on genomic variations involved in resistance 
development. However, data from high-incidence settings are only sparely available. Therefore we performed a 
systematic approach and analyzed a total of 97 M. tuberculosis strains from previously treated patients in Sierra 
Leone for mutations in katG, rpoB, rrs, rpsL, gidB, embB, pncA and where applicable in inhA and ahpC. Of the strains 
investigated 50 were either mono- or poly-resistant to isoniazid, rifampin, streptomycin, ethambutol and 
pyrazinamide or MDR and 47 fully susceptible strains served as controls. 

Results: The majority of isoniazid and rifampin resistant strains had mutations in tafG315 (71.9%) and rpoB531 
(50%). However, rpoB mutations in codons 511,516 and 533 were also detected in five rifampin susceptible strains. 
MIC determinations revealed low-level rifampin resistance for those strains. Thus, the sensitivity and specificity of 
sequencing of katG for detection of drug resistance were 86.7% and 100% and for sequencing of rpoB 100% and 
93.8%, respectively. 

Strikingly, none of the streptomycin resistant strains had mutations in rrs, but 47.5% harboured mutations in rpsl. 
Further changes were detected in gidB. Among ethambutol resistant strains 46.7% had mutations at emfc>B306. 
Pyrazinamide resistant strains displayed a variety of mutations throughout pncA. The specificities of sequencing of rpsl, 
embB and pncA for resistance detection were high (96-100%), whereas sensitivities were lower (48.8%, 73.3%, 70%). 

Conclusions: Our study reveals a good correlation between data from molecular and phenotypic resistance testing in 
this high-incidence setting. However, the fact that particular mutations in rpoB are not linked to high-level resistance is 
challenging and demonstrates that careful interpretation of molecular resistance assays is mandatory. In addition, 
certain variations, especially in gidB, appear to be phylogenetically informative polymorphisms rather than markers for 
drug resistance. 



Background 

With more than 9 million new tuberculosis (TB) cases 
and about 1.7 million deaths in 2009 [1] TB remains one 
of the most serious infectious diseases worldwide. Treat- 
ment and control of TB is further complicated by the 
emergence of drug resistant and even multi drug resist- 
ant (MDR) strains [resistance to at least isoniazid (INH) 
and rifampin (RIF)] [2]. Among high-incidence settings, 
Sub-Saharan Africa is eminently affected with two 



* Correspondence: sfeuerriegel@fz-borstel.de 

'Molecular Mycobacteriology, Research Center Borstel, Borstel, Germany 
Full list of author information is available at the end of the article 

(J Bio Med Central 



million new TB cases per year [3]. This study focuses on 
Sierra Leone, a high burden country with an annual TB 
incidence rate of 574 per 100.000 people and an annual 
mortality rate of 149 per 100.000 people. Treatment 
options are further hampered by the fact that 23% 
among previously treated TB patients in Sierra Leone 
suffer from an MDR M. tuberculosis strain [4]. 

Rapid detection of resistance is the key task to ensure 
an effective treatment of patients and also to avoid fur- 
ther spread of resistant M. tuberculosis strains. Molecu- 
lar assays that detect the genetic variants that mediate 
resistance constitute a rapid alternative to conventional 
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drug susceptibility testing (DST) and may even be per- 
formed directly on clinical specimens without culture 
[5,6]. Therefore it is essential to elucidate the genetic 
basis of clinical resistance and to correlate phenotypic 
and molecular resistance data. 

Resistance to INH is predominantly mediated by one 
mutation in the katG gene at codon 315 which results in 
the complete or partial loss of catalase-peroxidase activ- 
ity [7]. Further mutations in the promoter regions of 
inhA [8] and ahpC [9,10] are associated with INH resist- 
ance. Mutations responsible for RIF resistance are pri- 
marily located in the so-called rifampin resistance 
determining region (RRDR; codon 507-533 according to 
E. coli numbering system) of the rpoB gene which 
encodes the beta subunit of the RNA polymerase [11]. 
Resistance to streptomycin (SM) is mediated by muta- 
tions in different genes. Polymorphisms in rrs and rpsL, 
coding for 16 S rRNA and the ribosomal protein S12, 
respectively, are mainly responsible for high-level resist- 
ance [12]. Recently, the gidB gene, which encodes a 
7-methylguanosine methyltransferase specific for 16 S 
rRNA, has additionally been associated with SM resist- 
ance [13]. Resistance to ethambutol (EMB) is primarily 
mediated by mutations in the embB gene, coding for an 
arabinosyltransferase participating in mycobacterial cell 
wall synthesis, with codon 306 being most frequently 
affected [14]. Furthermore, mutations in other parts of 
embB (e.g. codon 406) [15] and upstream of embA 
[15,16] and in embC [16,17] are also involved in EMB 
resistance. Resistance to pyrazinamide (PZA) is known 
to be mediated by mutations occurring throughout the 



pncA gene, encoding a pyrazinamidase [18]. Resistant 
strains lack pyrazinamidase activity which is essential 
for pro drug activation. 

Since the frequency and combination of resistance 
mutations differs depending on the geographical setting 
in which the specific isolate is found [19,20], it is im- 
portant to analyze Mycobacterium tuberculosis complex 
(MTBC) strains from different regions and to determine 
putative setting specific molecular markers. However, up 
to now data about the accuracy of molecular diagnostic 
methods in high-incidence settings, and especially in 
West Africa, is only sparely available. Therefore we car- 
ried out a population based study, involving MTBC 
strains from Sierra Leone, to determine the genetic basis 
of first line drug resistance and to compare results from 
molecular and conventional drug susceptibility testing. 

Methods 

Mycobacterial strains and growth conditions 

A total of 97 MTBC strains isolated from previously 
treated patients in Sierra Leone were included in this 
study. All smear positive cases registered for re- 
treatment (failure after at least 5 months, relapses or 
treatment after interruption) between March 2003 and 
June 2004 in the Western Area and Kenema districts in 
Sierra Leone were recruited. From the strains analyzed 
50 were resistant to at least one of the following drugs 
INH, RIF, SM, EMB and PZA and 47 strains were fully 
susceptible (see Figure 1). From the panel of strains ana- 
lyzed, 74 were M. tuberculosis and 23 were M. africa- 
num strains. Primary isolation and cultivation was done 
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Figure 1 Overview of the antibiotic resistance profiles of the strains analyzed. A total of 97 M. tuberculosis and M. africanum strains from 
smear positive, previously treated patients from Sierra Leone was included in this study. Samples were collected in 2003 and 2004 in the Western 
Area and Kenema districts. Of the strains analyzed 74 were M. tuberculosis and 23 were M. africanum strains. Abbreviations: INH, isoniazid; RIF, 
rifampin; SM, streptomycin; EMB, ethambutol; PZA, pyrazinamide; R, resistance. 
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at the Supranational Reference Laboratory in Borstel as 
described previously [21]. 

Drug susceptibility testing 

Drug susceptibility testing (DST) to first-line drugs INH 
(0.25 and 1.0 ug/ml), RIF (20.0 and 40.0 ug/ml), SM (4.0 
and 8.0 ug/ml) and EMB (1.0 and 2.0 ug/ml) was per- 
formed in Borstel by using the proportion method on 
Lowenstein-Jensen (LJ) medium. In case of insufficient 
growth on LJ media, DST was done by applying the 
modified proportion method in the BACTEC 460 TB 
system according to the manufacturer's instructions 
(Becton-Dickinson). Drug susceptibility testing to PZA 
(100 ug/ml) was performed by using the BACTEC™ Pyr- 
azinamide (PZA) Drug Kit in the BACTEC 460 TB sys- 
tem according to the manufacturer's instructions. 
Determination of minimal inhibitory concentrations 
(MICs) was done by applying the modified proportion 
method in the MGIT 960 TB system (test concentra- 
tions were 1.0, 0.5, 0.25, 0.125 and 0.063 ug/ml for RIF 
and SM, 100.0, 50.0, 25.0, 12.5 and 6.3 ug/ml for PZA). 

DNA Isolation, PCR and sequencing 

DNA was isolated as described elsewhere [22] and amp- 
lified using the primers and conditions listed in Add- 
itional file 1. The PCR products were sequenced using 
an ABI 3130a:/ Genetic Analyzer (Applied Biosystems, 
CA, US) and the ABI BigDye Terminator kit v.1.1 
according to the manufacturer's instructions. The se- 
quence data was analyzed using DNASTAR Lasergene 
version 8.0, with M. tuberculosis H37Rv DNA as refer- 
ence sequence. 

All strains were sequenced in the predominant resist- 
ance determining regions (RDR) of katG (codon 315), 
rpoB (codon 507-533 according to E. coli numbering), 
rrs (nt. 1401-1402), rpsL (complete gene), gidB 
(complete gene), embB (codon 306) and pncK (complete 
gene). Strains resistant to the respective drug but not 
carrying a mutation in the RDR were sequenced in the 
complete gene. Strains resistant to INH with no muta- 
tion in katG were also sequenced in the promoter 
regions of inhK and ahpC Similarly, we extended our 
sequencing effort to embC and embA for EMB resistant 
strains without any mutations in embB. 

Results 

In this study a total of 97 MTBC strains with known re- 
sistance patterns to first-line drugs from Sierra Leone 
(see Figure 1) were sequenced in genes previously 
described to be involved in resistance development. 

The population structure of the strains was analyzed 
in a previous study [21]. Briefly, the 74 M. tuberculosis 
and 23 M. africanum strains belonged to eleven different 
genotypes. The population diversity was high with two 



M. africanum lineages (West African I, n = 6; West Afri- 
can II, n = 17) and nine M. tuberculosis lineages (Haar- 
lem, n = 14; LAM, n = 15; EAI, n = 4; Beijing, n = 4; 
S-type, n = 4; X-type, n = 1; Cameroon, n = 4; Sierra 
Leone I, n = 7; Sierra Leone II, n = 10). To determine if 
certain mutations appear genotype specific, the occur- 
rence of identified polymorphisms was correlated with 
the genotype of the respective strain. However, all muta- 
tions detected by sequencing analysis were found inde- 
pendently from their phylogenetic background (data not 
shown). A detailed summary of the sequencing data is 
provided in Table 1 and in Figure 2 (a-e). 

Among the INH resistant strains 71.9% (23/32) carried 
a mutation in katG at codon 315. Out of these, 21 dis- 
played a mutation in katG only, while two strains 
showed mutations at /catG315 with additional mutations 
at codon 291 and codon 471, respectively. One strain 
each carried a mutation at codon 300, codon 302 and 
codon 329. Two resistant strains displayed a mutation at 
codon 463, which is a phylogenetic SNP [23] and was 
therefore excluded from further analysis. Four of the 
INH resistant strains had no mutation in katG. However, 
sequence analysis of the intergenic regions of inhK and 
ahpC revealed polymorphisms in those areas. Two 
strains carried a mutation in inhK at position -15 and 
one strain in ahpG at -57. All of the 65 INH susceptible 
strains lacked mutations in katG. Thus for detection of 
INH resistance, sequence analyses of katG had a sensi- 
tivity and specificity of 86.7% and 100%, in the strains 
analyzed. 

Among RIF resistant strains, 50% (8/16) carried a mu- 
tation in rpoB at codon 531. The second most frequent 
mutation was found at codon 526 (37.5%). One RIF re- 
sistant strain each showed a mutation at codon 481 and 
at codon 533, respectively. Out of 81 RIF susceptible 
strains 76 did not have any mutation in rpoB. The 
remaining five susceptible strains displayed mutations at 
codons 511 (n = 1), 516 (n = 3) and 533 (n = l), respect- 
ively. Sequence analysis and drug susceptibility testing 
has been repeated for those five strains, confirming 
results of the first analyses. Determination of MICs 
revealed low-level RIF resistance (0.25-1.0 ug/ml) for 
those strains (see Table 2). Given that the strains show- 
ing low-level RIF resistance are assessed as susceptible 
by using standard DST, sequence analyses of rpoB had a 
sensitivity and specificity of 100% and 93.8% for detec- 
tion of RIF resistance, in the strains analyzed. 

To investigate the genetic basis of SM resistance, all 
strains were first sequenced in the rrs gene. As none of 
the resistant strains displayed a mutation in this gene, 
sequence analysis of rpsL was performed. Among all SM 
resistant strains 27.5% (11/40) carried a mutation in rpsL 
at codon 43 and 20% (8/40) showed a polymorphism at 
codon 88. The remainder of the phenotypically resistant 
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Table 1 Mutations detected in all strains analyzed 


Gene Mutation 


Strains phenotypically 
susceptible to INH (n=) 


Strains phenotypically 
resistant to INH (n=) 


katG wild type 


65 


4 d) 


Ala291Thr (get/act); Ser31SThr (agc/acc) 




1 


Trp300Cys (tgg/tgc) 




1 


Thr302Arg (acg/agg) 




1 


Ser315Thr (acg/acc) 




21 


Asp329 (gac/ggac frameshift) 




1 


Arg463Leu (cgg/ctg) (2) 




2 


Gln471Arg (cag/egg); Ser315Thr (agc/acc) 




1 




Strains phenotypically 
susceptible to RIF (n=) 


Strains phenotypically 
resistant to RIF (n=) 


rpoB wild type 


76 




Thr481 Ala (acc/gcc) (3) 




1 


Leu51 1 Pro (ctg/ccg) (4) 


1 




Asp516Tyr (gac/tac) (4) 


3 




His526Tyr (cac/tac) 




3 


His526Arg (cac/cgc) 




3 


Ser531Leu (tcg/ttg) 




8 


Leu533Pro (ctg/ccg) <4) 


1 


1 




Strains phenotypically 
susceptible to SM (n=) 


Strains phenotypically 
resistant to SM (n=) 


rpsi wild type 


56 


21 


Lys43Arg (aag/agg) 




1 1 


Lys88Arg (aag/agg) 


1 


8 
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15 (10 with mutations 
in rpsi) 
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1 (Lys43Arg in rpsL) 


Leu16Arg (ett/egt); Ala200Glu (gcg/gag) 
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Leu16Arg (ett/egt); Ala205Ala (gea/geg) 




1 (Lys43Arg in rpsi) 
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1 




Gly92Asp (gaa/gac); Ala205Ala (gea/geg) 




3 (Lys43Arg in rpsi) 


Serl OOPhe (tct/ttt) 




1 


Vail 1 0Va (gtg/gtt); Ala205Ala (gea/geg) 


4 




Ala138Val (gcg/gtg) 




2 


Ala200Glu (gcg/gag) + Ala205Ala (gea/geg) 




1 


Ala205Ala (gea/geg) 


18 


4 (3 with mutations in rpsi) 
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Table 1 Mutations detected in all strains analyzed (Continued) 

Strains phenotypically Strains phenotypically 

susceptible to EMB (n=) resistant to EMB (n=) 

embB wild type 82 4 (5) 

Met306lle (atg/ata) 6 

Met306Val (atg/gtg) 1 

Trp332Arg (tgg/cgg) 1 

Gln497Lys (cag/aag) (3) 1 

Gln1002Arg (cag/cgc) 13) 2 

Strains phenotypically Strains phenotypically 

susceptible to PZA (n=) resistant to PZA (n=) 

pncA wild type 84 3 

-11 (a/g) 1 

Ile5 (atc/atcc frameshift) 1 

Thr47Ala (acc/gcc) 2 

Tyr64 (tat/tata frameshift) 1 
Lys96Glu (aag/gag) 1 

Gin 1 41 (cag/cag acggcgccag (insertion of 1 
10 bp -» frameshift) 

Ala146Glu (gcg/gag) 1 

Gly162Asp (ggt/gat) 1 

Leu172Pro (ctg/ccg) 1 

Sequence analysis and drug susceptibility testing has been repeated for all strains showing discrepant results. 
Abbreviations: INH, isoniazid; RIF, rifampin; SM, streptomycin; EMB, ethambutol; PZA, pyrazinamide. 



Those strains were additionally sequenced in the promoter region of inhk and ohpC. Two of them carried a mutation in inhA at position -15 (C — *T) and one in 
ahpC at -57 (C^T). 

(2) Arg463Leu is a phylogenetic SNP [23] and was excluded for further analysis. 

(3) These mutations were identified by sequence analysis of the complete gene. 

(4) These mutations confer low-level RIF resistance. 

(5) These strains were sequenced in the complete embCAB operon. One strain carried a mutation in embC [Val981Leu (gtg/ctg}]. 



strains (n = 21) did not carry a mutation in rpsL. Among 
all SM susceptible strains (n = 57), one had the codon 88 
mutation in rpsL as well (confirmed when retested). De- 
termination of SM MIC showed no elevated MIC for the 
respective strain compared to the H37Rv control (see 
Table 2). Taken together, these data resulted in a sensi- 
tivity and specificity of the DNA sequencing of rpsL for 
detection of SM resistance of 48.8% and 98.2%, 
respectively. 

Additionally all strains were sequenced in gidB, In 
this very polymorphic gene 16 different mutations have 
been found, which occurred alone or in combination 
(see Table 1). Noticeable is the high number of phylo- 
genetic polymorphisms. The Leul6Arg (ctt/cgt) muta- 
tion was exclusively found in strains of the LAM 
genotype (n = 12). All strains belonging to the WA1, 
WA2 and Beijing genotypes displayed the Ala205Ala 
(gca/gcg) mutation (n = 27) and in all EAI strains a 
combination of the ValllOVal (gtg/gtt) and Ala205Ala 
(gca/gcg) mutations was detected (n = 4). The role of 
mutations in gidB for resistance to SM needs to be fur- 
ther investigated. 



Among all EMB resistant isolates 46.7% (7/15) carried 
a mutation in embB at codon 306. One EMB resistant 
strain was found to have a mutation at codon 332, one 
at codon 497 and two strains carried a polymorphism at 
codon 1002. In four EMB resistant isolates no mutation 
in embB was detected. Sequence analyses of embC and 
embA revealed a mutation in embC [Val981Leu (gtg/ 
ctg)] in one strain. All EMB susceptible strains (n = 82) 
had a wild-type embB sequence. Thus for detection of 
EMB resistance, sequence analyses of embB had a sensi- 
tivity and specificity of 73.3% and 100.0%, in the strains 
analyzed. 

PZA resistant isolates showed a wide variety of 
changes, distributed throughout the entire length of the 
pncK gene, including its promoter. Single nucleotide 
polymorphisms (SNPs) occurred in one strain each at 
position -11 bp, at codons 146, 162 and 172. In addition, 
insertions of single nucleotides leading to open reading 
frameshifts were detected at codons 5 and 64; an inser- 
tion of 10 bp after codon 141 led to PZA resistance in 
one strain. In three resistant isolates no mutation in 
pncK was determined. Among all PZA susceptible strains 
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Figure 2 Overview of mutations detected in all strains analyzed in a) katG, b) rpoB, c) rpsL, d) embB and e) pnclK. The height of the 
columns represent the number of strains analyzed. On the x-axis the different mutated codons are shown. The white areas of the columns 
represent the fraction of susceptible strains, whereas the black areas correspond to the number of resistant strains. Abbreviations: WT, wild type; 
singletons, various codons that are affected in one strain only. 



(n = 87), 84 displayed the wild type sequence, whereas in 
three PZA susceptible strains mutations were detected at 
codon 47 (n = 2) and at codon 96 (n = 1), respectively. Se- 
quence analysis and drug susceptibility testing has been 
repeated for strains showing discrepant results, however 
leading to unaltered findings. Determination of PZA- 
MICs (see Table 2) revealed slightly elevated MICs for 
the strains carrying the mutation at codon 47 (25.0 ug/ 
ml) compared to the H37Rv control, but an unaltered 
MIC for the strain carrying the polymorphism at codon 
96. These data resulted in a sensitivity and specificity of 
the DNA sequencing of pncK for detection of PZA resist- 
ance of 70% and 96.6%, respectively. 

A summary of all sequencing, DST, MIC and genotyp- 
ing data is provided in Additional file 2. 



Discussion 

In this study we carried out an in depth investigation 
of molecular resistance mechanisms by correlating par- 
ticular genomic variants with phenotypic resistance in 
clinical isolates from a high-incidence setting in West 
Africa. For INH and RIF there is a close correlation 
between data from molecular and phenotypic resist- 
ance testing for resistance determination in the strains 
analyzed. Sensitivity and specificity of sequencing of 
katG for detection of drug resistance were 86.7% and 
100% and for sequencing of rpoB 100% and 93.8%, 
respectively. 

Overall, the correlation between molecular and 
phenotypic resistance testing for the determination of 
SM, EMB and PZA resistance was lower. Although 
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Table 2 Determination of minimal inhibitory 
concentrations (MICs) of potential low-level resistant 
strains (to RIF, SM, PZA) 



strain 


mutation 


RIF MIC [ug/ml] 


4518/03 


rpoB Asp516Tyr (gac/tac) 


0.5 


5472/03 


rpoB Leu533Pro (ctg/ccg) 


1.0 


10011/03 


rpoB Asp516Tyr (gac/tac) 


0.5 


3736/04 


rpoB Leu51 1 Pro (ctg/ccg) 


0.5 


6467/04 


rpoB Asp516Tyr (gac/tac) 


0.25 


H37Rv control 


wild type 


0.25 


strain 


mutation 


SM MIC [ug/ml] 


6463/04 


rpsi Lys88Arg (aag/agg) 


0.5 


H37Rv control 


wild type 


0.5 


strain 


mutation 


PZA MIC [ug/ml] 


4724/03 


pncA Thr47Ala (acc/gcc) 


25.0 


4730/03 


pncA Thr47Ala (acc/gcc) 


25.0 


6467/04 


pncA Lys96Glu (aag/gag) 


12.5 


H37Rv control 


wild type 


12.5 



specificities of sequencing of rpsL, embB and pncK were 
high (96-100%), sensitivities were lower (48-73%) due to 
so far unknown resistance mechanisms. 

However, while our results in principle support mo- 
lecular resistance testing, the finding that especially in 
rpoB and also in pncK particular mutations are not 
linked to high-level resistance is challenging and demon- 
strates that careful interpretation of molecular resistance 
assays is mandatory. Therefore, studies targeting new re- 
sistance mechanisms should include valid phenotypic re- 
sistance data and, to our opinion, a comprehensive 
database on genetic variations in resistance genes and 
the correlation with phenotypic resistance is necessary. 
Furthermore, the level of resistance mediated by particu- 
lar mutations and the clinical consequences need to be 
thoroughly investigated. In addition, especially variations 
in gidB appear to be phylogenetically restricted rather 
than being involved in drug resistance development. 

In our study the most frequent mutation among INH 
resistant strains has been detected in katG at codon 315. 
This SNP has been observed in numerous prior studies 
[24,25] and has clearly been correlated with INH resist- 
ance by loss of catalase activity. In two strains, in 
addition to variations at /cafG315, mutations at codon 
291 and 471 were detected. However neither mutation 
has been described in the literature before and the 
/<b£G315 mutation therefore represents the likely mech- 
anism for INH resistance in these strains. The mutation 
at codon 300 observed in one strain in our study has 
been previously reported by Richardson and co-workers 
[26], where loss of this mutation has resulted in rever- 
sion of INH resistance in a previously drug resistant 
strain. The mutation at codon 302 as well as the 



insertion at codon 329 has not been described previ- 
ously. Since they are restricted to INH resistant strains 
in our highly diverse MTBC collection, they represent 
potential new INH resistance mechanisms. Experimental 
evidence is required to validate this hypothesis. 

Of the four INH resistant strains that displayed a wild- 
type katG sequence, three had mutations in the pro- 
moter regions of inhk and ahpC The inhh mutation 
has previously been described in the literature [24] as 
being the most common variation in the inhh promoter 
region related to INH resistance. Mutations in ahpC 
have been found before, however to our knowledge not 
at this position. In one of the resistant strains no muta- 
tion was found in neither the complete katG gene nor in 
inhh or in ahpG. This result suggests a so far unknown 
resistance mechanism as being responsible for INH re- 
sistance of this strain. 

Mutations in rpoB at codons 526 and 531 occur most 
frequently in the RIF resistant strains analyzed. Those 
SNPs are located in the RRDR and are well known for 
mediating resistance [27,28]. The mutation at codon 
481, which only occurs in one RIF resistant isolate, has 
to our knowledge not been described previously. 

The mutations at codon 511 (Leu— > Pro), 516 (Asp— > 
Tyr) and 533 (Leu — > Pro) conferred low-level resistance 
in agreement with previous studies [29,30]. It has been 
shown that various substitutions in the same codon lead 
to different levels of resistance. For example mutations 
at codon 516 can confer either low- or high-level resist- 
ance depending on the amino acid change [30]. Further- 
more, the phenomenon of RIF low-level resistance has 
only recently been described in a work by Van Deun and 
colleagues [31], where mutations at codon 511, 516 and 
533 have been found in strains tested susceptible by the 
radiometric Bactec 460 TB and Bactec 960 MGIT meth- 
ods. Our data confirm the existence of low-level RIF re- 
sistance mediated by specific mutations in rpoB that is 
not detected by standard drug susceptibility testing 
methods. However, MIC values, especially for the muta- 
tions at codon 516 and 533, are even lower (0.5-1.0 ug/ 
ml) than have been described in the literature. This fact 
may be due to the presence of further mutations in the 
operon or in other regions of the genome. 

In a recent study [32] the therapeutic challenge of 
low-level RIF resistance has been addressed and may, 
according to the authors, be overcome by the application 
of higher RIF doses (20 mg/kg) in treatment regimens. 
However, the clinical relevance and interpretation of 
these data is still not fully understood and needs further 
investigation in animal treatment models or clinical 
trials. 

Despite these discordant findings, we found a good 
correlation between the results from molecular and 
phenotypic testing for INH and RIF, as has been 
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observed in another study [33]. In fact, the strains ana- 
lyzed in this study predominantly harbour well described 
mutations which allows for the application of standard 
sequencing protocols or commercial line probe assays. 

The analysis of SM resistance mechanisms revealed an 
interesting observation. None of the SM resistant strains 
carried a mutation in the rrs gene, although those muta- 
tions have been described as main resistance mechan- 
isms that confer high-level SM resistance [12]. Instead, 
the SM resistant strains in our study population carry 
mutations in rpsL at codon 43 or 88 or at various 
codons in the gidB gene. The two mutations in rpsL 
have been described previously to confer high-level SM 
resistance [28,34]. Polymorphisms in gidB were reported 
to confer a lower level of SM resistance [13]. However, 
due to a number of phylogenetic polymorphisms in gidB, 
cautious interpretation of sequencing data is mandatory. 
Leul6Arg (ctt/cgt) has been described previously as 
phylogenetic marker for the LAM genotype [35], which 
could be confirmed in this study. Additionally, a syn- 
onymous SNP at codon Ala205Ala (gca/gcg) was identi- 
fied as being specific for the WA1, WA2 and Beijing 
genotypes, as well as a combination of Ala205Ala (gca/ 
gcg) and ValllOVal (gtg/gtt) was determined as phylo- 
genetically specific for strains belonging to the EAI 
genotype. These mutations in gidB occurred both in SM 
susceptible and resistant strains, affirming their role as 
phylogentic SNPs rather than markers for SM resistance. 
Polymorphisms in gidB probably playing a role in SM 
resistance, as they occur exclusively in SM resistant 
strains and do not coincide with mutations in rpsL, were 
detected throughout the complete gene (codons 34, 65, 
71, 88, 91, 100, 138, 200). However, the actual import- 
ance of these SNPs for SM resistance needs to be inves- 
tigated in further studies. 

Reasons for the absence of rrs mutations in the strains 
analyzed and the shift to mutations in rpsL and gidB are 
mainly unclear, but are in line with previous studies 
reporting a disequilibrium in the distribution of resist- 
ance conferring mutations in different geographical areas 
or among strains of different genotypes [36-38]. Our 
findings confirm that the performance of molecular 
assays that only target particular mutations can be influ- 
enced by the differential prevalence of particular muta- 
tions in a given geographical area. Therefore, strain 
diversity needs to be considered and investigated before 
the new implementation of molecular assays in a study 
region. 

Among EMB resistant isolates, the most frequent mu- 
tation affected codon 306 (Met/lie) of the embB gene. 
This mutation has been described in various studies as 
the main mutation mediating resistance to EMB [14,39]. 
The mutation at codon 497 has also been previously 
described in clinical isolates [40]. Moreover, both 



mutations have been shown to confer resistance by 
transfer in a wild type genetic background using allelic 
exchange experiments [41]. However, the authors con- 
clude that single mutations only modestly increase re- 
sistance to EMB and additional so far unknown 
mutations are necessary to cause high-level resistance. 

The mutations at codon 332 and 1002 determined 
here have not been described before. The impact of 
these changes has to be investigated in further studies. 
In four resistant strains no mutations were detected in 
the embB region analyzed. Additional analysis of the 
embC and embA genes revealed a mutation in embC 
[Val981Leu (gtg/ctg)] in one strain, belonging to the 
Haarlem genotype. However, this mutation has been 
described earlier as being specific for the Haarlem geno- 
type and is not associated with resistance to EMB [16]. 
As mentioned above, other so far unknown resistance 
mediating mechanisms are probably responsible for the 
resistance phenotype in these four strains. 

Mutations or insertions in the pncK gene are known 
to mediate PZA resistance [42,43], as observed in our 
study. No hotspot region has been determined, since 
polymorphisms occur throughout the complete gene. 
However, according to our data some specific mutations 
do obviously not mediate resistance that is detectable by 
applying standard critical concentrations. In the panel of 
strains analyzed, two susceptible strains carry a SNP at 
codon 47 and one displays a mutation at codon 96. 
PZA-MIC determination for these strains revealed 
slightly elevated values for the strains carrying the muta- 
tion at codon 47 (25.0 ug/ml) compared to the H37Rv 
control. In a recent study it has been shown that the site 
of the mutation is leading to varying efficiencies of the 
mutated pyrazinamidase mediating a wide range of re- 
sistance levels from low to high [44]. As the mutation at 
codon 47 has previously been described by Jureen and 
co-workers [42] in PZA resistant strains, further investi- 
gations are necessary to determine if additional muta- 
tions in other parts of the genome might be responsible 
for the observed low-level resistance in the strains ana- 
lyzed in this study. Out of all PZA resistant strains three 
carried the pncK wild type sequence. This indicates that 
further mutations in as yet unidentified genes are also 
important for mediating PZA resistance. 

Conclusions 

Although resistance mechanisms to INH and RIF are 
well understood, unknown resistance determining 
regions and resistance mediating mechanisms appear to 
play an important role for SM, EMB and PZA, where we 
observed a relatively low sensitivity for detection of re- 
sistance by analysis of common genes. Therefore, it is es- 
sential to gather information on further mechanisms 
leading to drug resistant MTBC strains. For the design 
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and implementation of molecular resistance assays it is 
fundamental to consider strain diversity with respect to 
resistance mutations in a given geographical setting. Fi- 
nally, it should be noted that not all variations in well 
described resistance genes are related to the development 
of high-level resistance, a finding arguing for a very care- 
ful interpretation of molecular resistance assays. 
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Additional file 1 PCR primers and conditions used for amplification and 
sequencing. In this table a summary of all primers, including 
oligonucleotide sequences, used in this study for both DMA amplification 
and sequencing is given. 

Additional file 2 Summary of all sequencing, DST, MIC and genotyping 
data. This table summarizes all data generated in this study. It comprises 
sequencing, DST (drug susceptibility testing) and MIC (minimal inhibitory 
concentration) testing results as well as all genotyping data. 
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